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Abstract 

We consider the phase-locked dynamics of a Josephson junction driven by finite-spectral- linewidth 
ac current. By means of a transformation, the effect of frequency fluctuations is reduced to an effec- 
tive additive noise, the corresponding (large) dephasing time being determined, in the logarithmic 
approximation, by the Kramers' expression for the lifetime. For sufficiently small values of the 
drive's amplitude, direct numerical simulations show agreement of the dependence of the dephasing 
activation energy on the ac-drive's spectral linewidth and amplitude with analytical predictions. 
Solving the corresponding Fokker-Planck equation analytically, we find a universal dependence of 
a critical value of the effective phase-diffusion parameter on the drive's amplitude at a point of a 
sharp transition from the phase-locked state to an unlocked one. However, for large values of the 
drive amplitude, saturation and subsequent decrease of the activation energy are revealed by simu- 
lations, which cannot be accounted for by the perturbative analysis. The same new effect is found 
for a previously studied case of ac-driven Josephson junctions with intrinsic thermal noise. The 
predicted effects are relevant to applications to voltage standards, as they determine the stability 
of the Josephson phase-locked state. 
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1 Introduction 



A particle in a viscous medium, in the presence of a spatially periodic potential, can be driven by 
a time-periodic force at a nonzero average velocity vq determined by the resonance condition, 

2'Km/ujQ^l/vo, (1) 

where loq is the driving frequency, an integer m is the order of the resonance, and I is the period of 
the potential. This phenomenon was experimentally observed in the form of phase-rotating states 
{Shapiro steps) in small Josephson Junctions (J J) driven by ac bias current [Q. Later, it was 
shown that the same effect can also be realized for a fluxon (magnetic flux quantum) moving in a 
periodically modulated long J J under the action of the ac bias current The latter effect has 
been recently experimentally observed, in the form corresponding to both the fundamental resonance 
[to = 1 in Eq. (^] and higher-order ones, in a long circular JJ with an effective spatial modulation 
induced by uniform dc magnetic field j^]. 

In terms of the phase difference 6 of the superconducting wave function across the junction, an 
ac-driven small J J is described by a pendulum equation (for a detailed derivation see the book Q) 

({^6 d9 

+ sin 61 = -a— + e cos {u^t + tpit)) + j{t) . (2) 

Here, time is measured in units of the inverse Josephson plasma frequency, a is the normalized 
JJ conductance, e is the ac-current amplitude normalized to the Josephson critical current, ujq is the 
normalized driving frequency, ijj is an arbitrary phase, and j{t) represents intrinsic thermal noise in 
J J. The equation of motion for the fluxon in the above-mentioned long circular J J in the magnetic 
field takes exactly the same form in the "nonrelativistic" limit, i.e., if the fluxon's velocity is much 
smaller than the Swihart velocity of the junction. The resonant relation (|l|), where I — 2ii, implies 
that, neglecting the noise, the ac drive in Eq. (^ may support phase rotation of the pendulum at 
the average velocity dO/dt = vq in the presence of the friction. Below, we assume that the spatial 
modulation period is always normalized so that I = 27r, i.e., vq = w^jm. We note that dO/dt is 
proportional to the voltage across J J, hence the phase rotation in the ac-driven J J gives rise to a 
nonzero dc voltage, a feature that is used in Josephson voltage standards (see, e.g., Ref. Q and 
references therein). 

In real applications, the driving ac signal is always slightly nonmonochromatic, having a finite 
width 5u} in the spectral domain; in other words, ?/; in Eq. (||) is not a constant phase, but rather a 
slowly varying function of time, representing random phase fluctuations of the driving signal. This 
introduces a finite lifetime H of the phase-locked state, which is of direct relevance to applications, 
affecting the stability of the Josephson voltage standards. 

Intrinsic thermal fluctuations, represented by the term j(t) in Eq. (^, also contribute to dephas- 
ing of the ac-driven motion. In terms of the small ac-driven JJ, the effect of thermal fluctuations was 
considered in earlier works where Eq. (H) with the monochromatic drive and thermal noise 

was reduced to a Langevin equation for a particle driven by a random force in a periodic potential. 
The phase-locked state is then represented by the particle trapped at a minimum of the potential, 
and the dephasing implies that the particle is extracted by the random force from the trapped state. 
The corresponding dephasing time was taken as the inverse Kramers' escape rate i.e., 

H ~ exp (AC//T) , (3) 

where AJ7 is the difference between maximum and minimum values of the effective potential, and 
T is the temperature (see exact definitions below). 

In this work, we focus on effects of the frequency fluctuations in the ac drive. In section 2, 
we demonstrate that Eq. (H) with a finite linewidth of the ac signal can be transformed, at the 
first order of the perturbation theory, into an equation driven by a strictly monochromatic signal 
and an additive random force. However, in contrast to the random force representing the intrinsic 
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thermal noise [j{t) in Eq.(g)], the correlator of the effective additive noise generated by the ac-drive's 
nonmonochromaticity does not contain the friction coefficient a, as this correlator, which has a non- 
thermal origin, does not obey the fluctuation-dissipation theorem. Next, using the energy-balance 
technique , in section 3 we reduce the monochromatically driven equation with the random force 
to a Langevin equation in a periodic potential U{6), which makes it possible to predict the dephasing 
time by means of the Kramers' expression (|3|). 

In section 4, we specially consider the situation when the small amplitude e of the driving signal is 
close to the phase- locking threshold ethr |§1 ■ By means of the Fokker-Planck equation corresponding 
to the above-mentioned effective Langevin equation ^ , we demonstrate that a critical value Slo 
of the driving-signal's linewidth, at which a sharp transition from the locking to unlocking (in the 
form of random 2tt phase slips) takes place, may be represented as a universal function of the drive's 
amplitude e. 

In section 5, we present results of direct numerical simulations of Eq. (|^) with the nonmonochro- 
matic drive, which are reported in the form showing the logarithm of the dephasing time H as a 
function of e and of the linewidth Slo. If e is above the threshold value ethr, and is not too large, 
the numerically found lifetime H is found to be quite close to that predicted by the perturbation 
theory. However, at large values of e the simulations reveal a new effect, which cannot be pre- 
dicted by the perturbative analysis: H(e) reaches a maximum value and then decreases. As the 
nonmonotonic character of the dependence H(e) and the existence of the maximum in it are quite 
important features, in section 6 we report results of direct simulations of the mode with the strictly 
monochromatic drive and intrinsic thermal (additive) noise. We conclude that the dependence H(e) 
in this case has the same nonmonotonic character. Although the latter model was studied earlier 
[|[ 0, this feature was not reported. 

It is relevant to mention that, besides small and long JJs, essentially the same dynamical model 
as the one considered in this work applies to ensembles of oscillators coupled via a mean field, which 
may be laser arrays or biological oscillators [pj| . As is known, the global coupling may synchronize 
the phase-rotation states of the oscillators, each of them being driven by the mean field. On the 
other hand, various perturbations affecting the mean field make it a slightly nonmonochromatic drive 
[HI |l3|. Thus, desynchronization of the globally coupled rotating oscillators is another manifestation 
of the problem considered in this work. 



2 Transformation of frequency fluctuations into an additive 
noise 

We begin our analysis from Eq. (^ without intrinsic thermal noise, i.e., with j = 0. White- noise 
fluctuations of the ac-drive's frequency, ujlt) = dTp/dt+coo, are assumed to be subject to the Gaussian 
correlations: 

Q being intensity of the fluctuations. The relative linewidth of the ac drive, Slo = {lu — luq) /wq 
(calculated at —3 dB level), can be then estimated as Suj w 2Ail/ujQ. We note that even low-quality 
sources of radiofrequency radiation, that may be used as the ac drive for JJ, have Suj ^ 10""^, while 
the dissipative constant a in JJ, although being small, is normally in the range a ^ 10~^, therefore 
we hereafter assume Suj ^ a. In other words, we may assume that a characteristic time of the 
variation of the random ac-drive's phase shift ^p(t) is much larger than the relaxation time l/a. 

To convert the frequency fluctuations into an effective additive noise, which is more convenient 
for the subsequent analysis, we define a new time variable which includes a slowly varying stochastic 
term, 

T = t + x{t),^ithxit) = ujo^m- (5) 
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Then, transforming the time derivatives d/dt into d/dr according to this definition, making use of 
the above relations ^cj <C a <C 1, and keeping newly appearing small perturbations at two lowest 
orders, we cast the underlying Eq. (||) into a form 

dO , , dx d^e dx d9 

_+sm0 = -a-+.cos(a;or)-2-. — -a-.-. (6) 

In fact, the last term in Eq. (||) is much smaller than the previous one, as a is small, and, in the 
first approximation, one may substitute, in the latter term, d^O/dr'^ « — s'mO. Thus, the final form 
of the perturbed equation, in which the frequency fluctuations were converted into the effective 
additive random force, is 



d'^e . ^ dO , , dx 

—r + smt' = -a— + ecos (luot) + 2 — 
dr'' dr dr 



2 . ^^..^ a— + ecos (cjqt) + 2— • sin0 . (7) 



The time transformation (|^) affects the Gaussian correlator (^. It is easy to find that, in terms 
of the renormalized time and renormalized random phase x (see Eq. (^)), an exact form of the 
correlator is 

/dx{r) dx{r)\-' dx{r') / dx{t' 



2n 6{t - 






1 - 


dr 



\ dr \ dr J dr' \ dr' 

However, in view of the smallness of the frequency fluctuations, in the lowest approximation we may 
adopt a simple form of the correlator. 



/dxir) dxiT')\ 2n 



(8) 



Equations (0) and (|^) will be a basis for further analysis, while numerical simulations will be run 
for the full underlying equation (0). 



3 An effective Langevin equation and estimate for the life- 
time of the phase-locked state 

In the zero-order approximation, which implies e = a = dx/dr = 0, Eq. has a known solution, 

(?oW = 2am((r-ro)/fc;fc), (9) 

corresponding to the phase rotation at a nonzero average frequency (phase velocity) 

wo = T:/kK{k) . (10) 

Here, am is the Jacobi's elliptic amplitude with the modulus k {0 < k < 1 ), K{k) and E{k) 
are complete elliptic integral of the first and second kinds, and tq is an arbitrary constant. In 
this approximation, Eq. (||) conserves the energy, the value of which for the law of motion (||) is 
determined by the value of fc, 

A possibility to support persistent phase rotation at a velocity ujq in the presence of the friction 
by the monochromatic ac drive is predicted by the energy-balance equation pi. To this end, we 
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calculate the net rate of the change of energy due to the action of the friction and drive, averaged 
over the rotation period 2it/ujo, under the resonance condition given by Eq. (|l|): 



(12) 



where the overbar stands for the time average, tq being the same constants as in Eq. (^), and uji is 
the amplitude of the resonant harmonic in the Fourier decomposition of the time-dependent velocity 
(instantaneous frequency) dOo/dr, taken as per the unperturbed law of motion (^. An elementary 
calculation yields 



d9o 
dr 



UJl 



2Q 



1 + Q2 



Q = exp 



^ttK 



(v^ 



/K 



(13) 



~ k^K{k) ' 

{Q is called the Jacobi parameter) 

dS/dr — 0, is possible at two constant values (one stable and one unstable) of the phase difference 
between the ac drive and the rotating pendulum. 



Then, the phase-locked ac-driven regime, corresponding to 



Wo To 



± cos ^ 



eujQLJi \ dr 



(14) 



[||, provided that the amplitude e exceeds the threshold value 



ethr 



uJoLOi \ dr 



(15) 



In the presence of the additive random force in Eq. (j^), a perturbed equation of motion for 
the rotating pendulum can again be obtained from the averaged energy-balance equation, which 
has the same form of Eq. (12) with a difference that now to is a slowly varying function of the 
time T [roughly speaking, varying as slowly as the random phase ?/;(<:) in the underlying equation 
(^]. Notice that the time dependence of To determines a change Sujq = wo {dTo/dr) of the average 
phase-rotation velocity, and the latter may be related to the change of the energy (|l^), through 
its kinetic part, as d£ — S'Sluq, where £' stands for d£ /dujo, calculated for the unperturbed law of 
motion given by Eq. (|), £' = 4,K'^ l[Kk{K + kdK/dk)]). Thus, dsjdr^ S'ujq (dVo/dT^), and the 
balance equation (|l^) takes the form 



LUqS 



,cPto 

dT^ 



eujQUJi cos (wqTo) 



2u;o^-sin(0(T)), (16) 



the second term in {l + 2dTo/dT) being a contribution to the energy dissipation rate due to the small 
change of the average velocity, while a similar correction to the last term in Eq. (16) is negligible. 
Eq. dig) can be transformed into a more convenient form by defining ujqTo{t) = C(t) : 



dr"^ 



2a 



Mo 
dr 



dC 
dr 



euf^uji cos C] -f 2ti)Q sin {9{t)) 



dx 
dr' 



(17) 



Thus, we have arrived at an effective Langevin equation Q for a particle driven in a viscous 
medium by the sum of a regular force, represented by the terms in the square brackets, and a 
stochastic force, represented by the last term in the equation. The subsequent step is, following a 
well-known procedure ||l^, ^, to introduce the Fokker-Planck (FP) equation corresponding to this 
Langevin equation, taking into regard the correlator (||) . An essential feature of the thus derived FP 
equation is the presence of the extra multiplier sin^ [^(t)] in front of the diffusion (second-derivative) 
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term in it, due to the multiplier sin [9{t)] in the stochastic-force term in Eq. (|l7|). The coefhcient 
sin^ [9{t)] may be averaged in time as per the unperturbed law of motion @. It is easy to calculate 
the average value: 



sm 



2-P 



m 

K{k) 



(18) 



[note that, in the limit ^ 0, the expression (|l8|) approaches an obvious value 1/2]. 

The FP equation takes essentially the same form as it would take in the known problem [|[ Q 
of the dephasing of the ac-driven JJ phase rotation under the action of the additive thermal noise 
represented by the term j{t) in Eq. (^) [the most important difference of the effective Langevin 
equation ( [l7| ) from its counterpart in the thermal-noise problem is the presence of the multiplier 

With the FP equation taking the usual form, one 
to determine the lifetime of the ac-driven state, as 
it was done for the case of the thermal noise in Refs. 0]- In particular, an effective potential 
corresponding to the potential force, i.e., the expression in square brackets in Eq. (|l7|), is 



sin [0{t)] in the last term of the equation 
can directly use the Kramers' expression ( 



U (C) = —eugOJi sin C -I- uiQa 



dr 



c, 



hence the potential-barrier height At/, which should be substituted into Eq. (^, can be easily 
found as a difference of values of the potential taken between two points where the above-mentioned 
potential force vanishes. As a result, we obtain 



\/l - (ethr/e)^ - (ethr/e) cos ^ (ethr/e) 



(19) 



where the definition (g_5|) for the threshold value of the amplitude has been used to simplify the 
expression. 

However, an additional difference of the present case from the thermal-noise problem, that must 
be taken into regard before applying the expression (||), is that the frequency fluctuation intensity 
appearing in the correlator (|^) does not obey the fluctuation-dissipation theorem, and therefore it 
does not include the intrinsic dissipative constant a of the pendulum (JJ). By properly defining the 

effective temperature Toff = 4wQsin^ 9q{t) fl/aes, where acs = 2a(d0o/dT)^ is the effective friction 
constant from Eq. (|l7|), and using the potential-barrier height AU (|9[) and the expression (|8|), we 
can rewrite the Kramers' expression (^) as: 



exp 



Qeak'^EQ 



[(2-A:2)£;_2(l-fc2)X](l 



1 - (ethr/e) - (ethr/e)cos ^ (ethr/e) 



(20) 



This is the main prediction of the analytical consideration, which will be compared to results of 
direct simulations of Eq. (^) in section 5. 



4 Dephasing the phase-locked state near the locking thresh- 
old 

In this section we investigate the system described by the Langevin equation ([l7|), by explicitly solv- 
ing the associated FP equation. However, we can first simplify Eq. (17), recalling the fundamental 
physical condition according to which the frequency fluctuations are much smaller than a, or, in 
other words, the random force varies on a time scale 3> 1/a . Consequently, the acceleration term 
on the left-hand side of Eq. ( p^ ) may be neglected as compared to the velocity term, which yields 
a simplified Langevin equation 
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C =-Fo + FiCosC + f{t), 



(21) 



where Fq = uJo/2, Fi = eujQU!i[2a{d9o / dt)'^] ^, and 



fit) 



d0o 

dr 



^•sin[0(r)] 



(22) 



The FP equation (in this case, it is, in fact, the Smohichowski equation |T^) for a probabihty 
distribution function P(C,i), corresponding to Eq. ( pl| ) with the Gaussian correlator (||), is 

Pt=FismC-P+ {Fo - Fi cos C) Pq + SP^c , (23) 
where the subscripts stand for the partial derivatives, and 



S= sin^0o(T) 



dr 



(24) 



where Q. is the same as in Eqs. (^) and (||), and the average value sin^ Sq{t) is given by Eq. 

Information about the distribution of the phase C can be obtained from the stationary version 
OfEq. dH), 

d"^ P dP 
D—^-{\-bcosO—-bs\nC-P, (25) 

where the final set of notation is b = Fi/Fq and D = S/Fq. These two parameters are interpreted 
, respectively, as the ratio of the drive's amplitude to the friction coefficient, and an FP diffusion 
coefficient, which is proportional to an effective drive's linewidth. 
In the absence of the diffusion (Z? = 0), a solution to Eq. (E5h is 



P(C) = (27r)"^ Vl -^V(l - &cosC), 



(26) 



where the normalization J^^ P{(^)d( = 1 is imposed. The solution (|2^) is regular at 6 < I, while 
the singularity at 6 = 1 exactly corresponds to the drive's amplitude attaining the threshold value 
(p^, i.e., to the onset of the phase-locking regime. 

Collecting results produced by the numerical solution of Eq. (25), we have concluded that it is 
possible to define a critical value Da oi D, at which a sharp transition from the phase-locked state 
at D < Dcr to an unlocked one at D > Da takes place. The transition is still better illustrated 
by consideration of the probability fiux, J = — [DP' + (1 — fecosC)P], in terms of which the time- 
dependent FP (Smoluchowski) equation ( p3| ) is written as Fo^^Pj-l- dJ/dC, — 0. At the points of 
a minimum of the stationary distribution function, \ J\ gives the phase-slippage rate, i.e., a rate of 
the transition from a vicinity of a phase-locked point to a point differing by a phase shift ±27r. In 
Fig. 1, we display the phase-slippage rate vs. D at different constant values of b, as found from the 
numerical solution of Eq. (p5|). The existence of critical values _Dcr, such that virtually no phase 
slippage takes place at D < Dcr^ is evident. Note that for & < 1, when the locking is impossible even 
to the monochromatic drive in the absence of the noise {D = 0), Dcr does not exist. For 6=1, i.e., 
exactly at the locking threshold (|l^ ), Dcr = 0, which cannot be seen on the logarithmic scale used 
in Fig. 1. It is necessary to mention that a picture which may be interpreted as showing the fiux J 
as a function of 6 at different fixed values of D is available in the book nevertheless, we find it 
relevant to present Fig. 1 here, as we need to display J{D) at different fixed values of b (otherwise 
the existence of the critical values Z?cr is not obvious). 

The sharp unlocking transition can also be seen in terms of the ratio of the aforementioned 
minimum value of P(C) to its maximum value, which is attained fairly close to the unperturbed 
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(i.e., pertaining to the monochromatic drive) locking point. These data (not displayed here) show 
that the ratio is virtually equal to zero at I? < Dcr, and abruptly begins to increase exactly at 
D = Dci- In order to quantify Dcr we define it as the value of D for which J ~ 10^^. 

Figure 2 shows the most important characteristic of the unlocking transition, viz., Dcr vs. the 
effective drive's amplitude b. The dependence is nearly linear, except for the region < 6 — 1 <C 1, 
i.e., just above the locking threshold (|l5|). It is difficult to show the dependence in this region 
directly, therefore, in Fig. 2 we instead display r{b) for 6 — 1 <C 1, where r is defined so that Da is 
approximated by an expression Ci(6 —1)'' with a suitable constant Ci. 

An asymptotic value of r for 6—1-^0 can be found analytically, as one can use the exact 
solution (|2^) to describe an approximate form of the distribution function, except in a sensitive 
region of small C: P{C) ^ l/sin^(C/2). At small one can expand Eq. (|25|), taking into regard 
that b — 1 and D are now small too. This yields an equation, 

d^P /I 2 b-l\ dP 

where r] = D~^/^C,. Thus, the solution in the sensitive region (which is 77 ^ 1, or C ^ \/6 — 1 ^ D^^'^) 
depends on the single parameter (1 — 6)/_D^/^. Although a solution to Eq. ( p7| ) can be matched 
to the aforementioned approximation P(C) ^ l/sin^(C/2), valid at larger C,^ only numerically, it is 
obvious that, as 6 — 1 ^ 0, the dependence Dci{h) must be — Ci{b ~ 1)'^/^, with Ci ~ 0.17 
found from numerical data. The value r — 3/2, obtained for 6 — 1 ^ 0, is in good agreement with 
the numerical results displayed in Fig. 2. 

The linearity of the dependence Dci {b) at large b can be explained in a very simple way: neglecting 
in this case 1 in the expression (1 — 6cosC) in Eq. (p5[), we immediately conclude that the asymptotic 
solution depends on the single parameter b/D, hence the dependence must take the form Dd- — C26, 
with a constant C2 ~ 0.14 found numerically. The latter result can also be interpreted in another 
way: the minimum (threshold) value of the ac-drive's amplitude necessary to support the rotation 
of the pendulum grows nearly linearly with the linewidth Slo, so that it may be approximated by 
ethr (Suj) ~ e[h|.(l + const-(5w), where e[|||, is given by Eq. ( p^ ) and the constant is roughly I/C2. 

These analytical results, which comply well with the numerical findings, justify the introduction 
of the very concept of the critical value Da- of the phase-diffusion constant in the Smoluchowski 
equation, which was originally defined above in a phenomenological way, just by looking at Fig. 1. 



5 Numerical results 

To check the limits of validity of the analytical results obtained above, we have performed numerical 
simulations of the full stochastic equation (|^), using a simple Euler scheme (the use of this scheme 
is considered in the book [^). As usual, we halved the time step until the results would converge 
to a stable value within few-per-cent accuracy. In this section, we will first focus on the case of 
the parametric noise, so we now set j{t) = 0, keeping the random term ip{t) in Eq. (H). The basic 
phenomenon sought in the previous analysis was the escape from the state synchronized with the 
external drive as per Eq. (^. However, since the J J is driven by the ac-term alone, once the system 
is no longer phase-locked, it cannot sustain progressive motion and will therefore quickly decay to 
the zero- voltage state. So the prediction of Eq.(^ ) actually refers to the lifetime of the phase- locked 
state; an abrupt transition from this state to the zero-voltage one (after about 300 time unit) is 
evident in Fig. 3, which displays the time dependence of the phase velocity, found in a typical run 
of the simulations of the stochastic equation (^ . 

To estimate the lifetime of the phase-locked state at a given "temperature" (spectral linewidth of 
the drive), we have run the simulations for many different realizations of the random phase ipit), and 
averaged the results for the lifetime. The number of the realizations was determined by the condition 
that the computed average has to converge to an established value. In Fig. 4, the logarithm of the 
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thus computed average lifetime is plotted versus the inverse linewidth, so that the potential barrier 
could be estimated from the slope of the linear part of this dependence. 
This method closely follows that of Ref. j^; the main difference is, as already mentioned, that 
we are not looking for mere phase-slippage, but for a jump to the state with zero average velocity. 
As it is clearly seen in Fig. 3, this occurs at a somewhat later time than the phase slips commence, 
although the difference is, typically, small (for instance, it is less than 20 time units in the example 
shown in Fig. 3). 

The numerically computed effective energy barrier (represented by lines with symbols) and the 
one predicted by Eq. ( p^ ) (the lines without symbols) are shown in Fig. 5(a). Taking into regard 
that no fitting parameter was employed, the agreement is very good near the threshold value (p^), 
which is ethr ~ 0.16 for the values of the parameters corresponding to Fig. 5(a). However, a drastic 
deviation from the analytical prediction is evident at larger values of the ac-drive's amplitude: while 
the analytical formula ( |20| ) predicts an almost linear increase of the effective barrier height with the 
ac drive amplitude, the numerical results show a maximum followed by a substantial decrease of 
barrier's height. 

It is relevant to mention that the nearly linear increase of the barrier height with e was predicted 
by the power-balance approach, which was employed above for the analytical consideration (see also 
Ref. 1^), a different method, based on an harmonic expansion, would result in a nonmonotonic 
dependence of the barrier height on e ^. In the case of the additive noise, the latter method 
produced the energy barrier demonstrating a Bessel-functional behavior, typical of the rf-induced 
current steps in JJs In our case, however, such a dependence cannot be analytically justified. 

Moreover, even if the results of Ref. ^ are formally applied to our case, yielding 



where Ji is the Bessel function, the maximum of AU would occur at a much larger value of 6 [6 ~ 45, 
instead of 4 in Fig. 5(a) for a = 0.01]. Thus, the phenomenon reported here is an essentially new 
one, and still remains to be explained. 

It should, moreover, be noticed that the Bessel-function approximation similar to Eq. ( |2^ ) gives, 
according to Ref. ||^, a good estimate for the energy barrier only in the limit << 1, deviations 
occurring already for uJq'^ ~ 0.05. Therefore, it is not surprising that, for the parameters considered 
here {lOq'^ = 0.25), the agreement with the Bessel-function behavior is very poor. 

To check that the newly found dependence of the effective barrier height vs. the normalized 
drive's amplitude is not due to some particular feature of the parametric noise, we have also per- 
formed extensive simulations of the same stochastic equation (^, but with the additive noise and 
strictly monochromatic driving signal, for the same values of parameters as those used above in 
the case of the frequency fluctuations. A typical example is shown in Fig. 5(b), together with the 
theoretical estimate of the energy barrier according to Ref. |Q, ^ : 



In this case too, a strong deviation of AU (b) from the linear dependence occurs at relatively low 
values of the drive amplitude (& ~ 7), although they are higher than those in the frequency-noise 
case [which are 6 ~ 4, see Eq. 5(a)]. 

6 Conclusions 

Results reported in this work are of relevance for applications to systems in which Josephson junc- 
tions are phase-locked to an external ac source, such as voltage standards. By substituting reasonable 
experimental values into Eq. (|20|), we can estimate that an ac source with the relative linewidth 
better than 10~^ is needed if the lifetime of the order of 1 second is required for the measurement 
system. 




(28) 




(29) 
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We also note that our approach could be applied to another problem: a pendulum driven by 
an ac signal whose frequency is subject to a systematic (rather than random) change, i.e., a zero- 
linewidth but variable-frequency drive. A system of the latter type was considered in Ref. for 
a soliton in a perturbed nonlinear Schrodinger equation. 

This work was performed in the framework of the bilateral-cooperation agreement between Con- 
siglio Nazionale di Ricerca (Italy) and the Israeli Ministry of Science and Technology, under the 
project Nonlinear Dynamics of Josephson Networks. The authors appreciate useful discussions with 
G. Costabile and access to computer facilities at the University of I'Aquila (Italy). 
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Figure Captions 



Fig. 1. The phase-slippage rate J vs. the diffusion parameter D at different fixed values of the 
normalized drive's amplitude b, which are indicated near each curve. 

Fig. 2. The normalized noise threshold Dcr vs. the normalized drive's amplitude b (solid line). The 
dashed line shows the dependence of the exponent r defined in the text. 

Fig. 3. A typical example of the evolution of the phase velocity d{t) obtained from the numerical 
integration of the stochastic equation (H). A loss of phase-locking occurs at t ~ 300. The parameters 
are a ^ 0.01, e = 0.2, n = 0.01, ujq 2. 

Fig. 4. Lifetime of the phase-locked state, plotted on a logarithmic scale versus the inverse linewidth. 
The parameters are a — 0.01, wo — 2. The estimated slope is 0.00031 for e = 0.2, and 0.0012 for 
e = 0.5. 



Fig. 5. Dependence of the normalized energy barrier AU on the normalized drive's amplitude 6, 
for the case of frequency fluctuations (a) and additive noise (b). Symbols represent the energy 
barrier estimated numerically on the basis of Eq. (^) for ujq — 2 and a — 0.1 (squares) and 
a = 0.01 (triangles), connecting lines being a guide for the eye. The curves without symbols 
represent analytical predictions, viz., Eq. ( p0| ) for the frequency-fluctuation case [in the panel (a), 
the continuous and dashed curves pertain, respectively, to a = 0.1 and a = 0.01], and Eq. (29) for 
the additive-noise case. 
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